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I. INTRODUCTION 



In supersymmetric (SUSY) models without dimension-3 supersymmetry-breaking op- 
erators, gauginos are massless at the tree level and obtain non-zero masses solely from 
radiative corrections Jl],|2|||. This means that the gluino is light and the lightest neu- 
tralino is nearly a pure photino. Farrar {|,|5]|| found that the light gluinos and photinos 
arising from this scenario are consistent with the present experimental constraints.^ Al- 
though it was once generally believed that the light gaugino scenario conflicted with 
the cosmological relic abundance constraints, Farrar and Kolb showed that the pre- 
vious constraint calculations neglected the reaction channels which really control the 
relic abundance. Indeed, based on some simple estimates, they concluded that a light 
photino (the relic stable particle in the present SUSY scenario) might be a significant 
dark matter candidate. However, their estimates were based on the approximation that 
only a single reaction dominates the relic abundance evolution and that the abundance 
evolution stops exactly when the dominant reaction rate becomes less than the Hubble 
expansion rate (the "sudden" approximation). 

In the present paper, we calculate the cosmological constraints for this scenario of 
light gluinos and photinos more carefully by integrating the Boltzmann equations for 
the relic abundance. The three most important reactions determining the photino abun- 
dance are R ^ «-> 7r 7, R° — > 7r + 7r~7, and R°R° — > X. The first two are related by 
crossing symmetry. In the limit that left and right handed squark masses are equal, such 
that charge conjugation is a good symmetry of the theory, and ignoring the momentum 

2 The recent ALPEH claim to exclude light gluinos JtJ assigns a la theoretical systematic error based 
on varying the renormalization scale over a small range. Taking a more generally accepted range of 
scale variation and accounting for the large sensitivity to hadronization model, the ALEPH systematic 
uncertainty is comparable to that of other experiments and does not exclude light gluinos [pf. 
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dependence of the matrix elements, this crossing relation indicates that both reactions 
may play an important role instead of one reaction dominating over the other. We use 
the results of our model calculations to help identify the cosmologically most promising 
values for phenomenologically important parameters such as the R° lifetime, which can 
help in laboratory searches. 

Let us now briefly introduce the relevant features of our SUSY scenario. Supersym- 
metric models with acceptable SUSY breaking phenomenology are generically invariant 
under a global chiral symmetry called i?-invariance. i?-invariance is broken sponta- 
neously by the vacuum expectation values of the Higgs fields associated with electroweak 
symmetry breaking, and by tree- level gaugino masses if they are present. i?-parity is the 
possible discrete remnant of this broken continuous symmetry. Under i?-parity, the 
gluino, photino, and squarks are odd, while ordinary particles (e.g., gauge and Higgs 
bosons and quarks) are even. .R-parity, which we shall assume is an unbroken symmetry, 
ensures that the lightest i?-odd particle is stable and prevents unacceptably rapid pro- 
ton decay. Thus, in calculating the relic density in SUSY, one first identifies the lightest 
.R-odd particle which usually is the lightest supersymmetric particle (LSP). Although 
the gluino may be the lightest particle in our scenario, it cannot exist in isolation today 
because it is not a color singlet. Bound to a gluon or a color-octet system of quarks 
and/or anti-quarks, it forms a color singlet hadron. The lightest of these is expected to 
be a gluon-gluino bound state called R°, whose mass should be comparable to that of 
the lightest glueball 0, [|. Because this is most likely heavier than the photino, it is 
the photino which acquires the role usually taken on by the LSP even though it may be 
heavier than the gluino.0 

3 In some SUSY-breaking models only the gluino is massless at tree level, while other gauginos have 
large masses. In this case the R° could be the LSP and relic R°'s would be the SUSY dark matter 
candidate. The dark matter density can be approximated as in [fta, accounting for only the R°R° self- 
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Tabic I: A list of SUSY mass parameters and ranges used in the analysis. 

Particle Mass Notation Min.(GeV) Max.(GeV) 




M S 



M 



rn 



0.2 



50 



1 



300 



1.4 
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Since freeze out occurs after the color confinement phase transition, only gluinos 
bound in color singlet states are relevant to our calculation. Among the bound states 
containing a gluino (i?-hadrons), the R° is expected to react most prominently with 
the photino because other i?-hadrons are significantly heavier and thus Boltzmann sup- 
pressed at the relevant temperatures. Furthermore, most of the other i?-odd states will 
contribute to the photino abundance only after having decayed to an R° channel. Thus, 
the photino relic abundance will be determined primarily by the reactions involving an 
R°, a 7, and non-SUSY particles. 

In our scenario, the photino abundance depends crucially on interactions of hadrons 
after the confinement phase transition, causing complications distinct from conventional 
scenarios where the freeze out occurs above the confinement transition temperature. In 
particular, we are only able to make reasonable guesses for the relevant reaction rates 
because of incalculable long-distance QCD effects and our lack of direct experimental 
data for the reaction rates of interest. Fortunately, we are still able to make useful 
predictions regarding the R° and 7 masses and R° lifetime. 

The relic abundance of photinos depends mainly upon the masses of 7 and R°, the 

cross sections for i?°7r — > 71-7 and R°R° — > X (where X denote any strongly interacting 

annihilation. This gives U, R oh 2 < 10~ 7 . That is, due to their strong interactions, R°'s stay in thermal 
equilibrium too long for their abundance to freeze out at a non-negligible value. Thus such SUSY- 
breaking scenarios do not provide a natural visible sector dark matter candidate unless the gravitino 
has acceptable properties. 
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light species of particles such as the pions), and the decay rates for R° — > ^nn and 
R° —> jit. The mass parameter space that will be explored in this calculation is justified 
in Refs. and || and is similar to that discussed in Ref. ||. The relevant mass 
parameters and their plausible ranges are shown in Table I. The gluino mass itself is 
unimportant, except insofar as it influences the R° mass, M. The relevant squark mass 



denoted Ms is a charge-weighted average of up- and down-squark masses. See Ref. JlO 
for squark mass limits in the light gluino scenario. 

In order to express some of the formulae showing numerical estimates concisely, we 
also define the following dimensionless ratios: 

m Ms M 



Z^S = 777TTTTF' 



0.8 GeV ^ 100 GeV m 
As pointed out in Ref. , the relic abundance is particularly sensitive to the parameter 
r. Note that because of the ranges we adopt, given in Table I, the range of r we explore 
is constrained for a fixed value of m. 

In the next section, we discuss the Boltzmann equation and some simplifying as- 
sumptions used to calculate the present photino abundance (density). In Section III, we 
briefly describe the reactions that are included in the simplified Boltzmann equations. 
The results of the integration are presented and analyzed in Section IV. In Section V we 
develop an effective Lagrangian description of the interaction between R°, 7 and pions 
which embodies the symmetries of the underlying theory as well as the crossing and 
chiral-perturbation theory constraints. Ignoring the possibility that a nearby R v reso- 
nance produces a strong momentum dependence, the two dominant reactions controlling 
the 7 abundance are determined by a single parameter. Using this approximation, we 
obtain an estimate of the cosmologically favored lifetime range of the R° as a function 
of its mass. We summarize our results in Section VI. In an appendix, we analyze the 
possible resonance enhancement of the i?°7r — > 777 cross section using a Breit-Wigner 



model. 



II. THE BOLTZMANN EQUATIONS 



The standard method of calculating the relic abundance is to integrate a simplified 

form of the Boltzmann equations 0,0]. We now briefly remind the reader of the general 

formulation. One can write the Boltzmann equations^ for the evolution of the particle 

density rij as 

dn, _ (W jA .._ B ..V n ^) 
—± + 3H nj = - y N l^f^i L 

x L T[n k Y[ n? - nf ft n? ft n x ) (2) 

where H is the usual Hubble expansion rate, Aji and Bji are sets of particle species 
relevant to the evolution of species j, and the summation is over all the reactions of the 
form jAji — > Bji. We have defined the thermal averaged transition rate as 



[[dp}(2nr6^(y Pk )\T jAj ^ Bj fexp(- £ E x /T) 

\ k J Aej.Aji 



ww""*^ c^s — (3 > 

Ti?^^/ ^ ex P(-^/ T ) (5) 

where n^ 9 is the equilibrium density^] g x counts the spin multiplicity, and IT^^^.J 2 

represents the spin averaged transition amplitude squared.^ In the case of one initial 

4 As usual, we have used the assumption of molecular chaos to obtain a closed set of equations. 
5 As usual, the particle described by this equilibrium density is assumed to have mass much greater 
than the temperature. 

6 The W symbol actually represents the number of transitions per unit time when all the initial state 



state particle, {W jA ji ^B ji V n ^ Aji ' > ) evaluates to a decay rate, whereas in the case when 
there are two initial state particles, (WjA ji ^B ji V n ( Aji ' > ) evaluates to the familiar (av) of a 
scattering reaction. For example, in the case of a photino density evolution determined 
only by the reaction 77 «-> X, the density labels become j = 7, Aji = A~ = {7}, 
and Bji = B~~ = {X}; the summation in % reduces to a sum over one element (the 
annihilation channel); and the transition rate per unit fluxes becomes (Wzz_+ x V) = 
{W(77 — ► X)). With the usual assumption that the final products X are in equilibrium, 
Eq. (H) reduces to the familiar equation (see, for example, pg. 120 of Ref. fll2"l ) 

^ + 3Hn~ = -(w(77 - *)>(n| - 2 ). (6) 

Note that Eq. (§) assumes that the fluid is rare enough to disregard degenerate pressure 
effects and assumes that time reversal is a good symmetry. More specifically, time 
reversal symmetry is encoded in the following identity used in obtaining Eq. @: 



Before we can utilize Eq. @ to determine the relic abundance of the photinos, we 
need to specify our model of H and the reactions that are involved. Because the universe 
is radiation dominated for the temperatures of interest, the equation of state is taken 
to be 3 x (pressure) = (energy density) and any possible spatial curvature is neglected. 
We also use equilibrium statistics with the number of relativistic degrees of freedom set 
to g* = 10.75. The resulting equation for the Hubble expansion rate as a function of 



temperature is H = y 87r 3 g*/90(T 2 /m p z), where m p i is the Planck mass. The reactions 
that can enter the Boltzmann equations include 7-R <-> X, 77 <-> X, R°R° <-> X, 
77r <-» R°, 77T7T <-> R°, and 771 i?°7r (X's denote any allowed light products that 
interact strongly or electromagnetically). 

reactants have equal unit flux. The V symbol represents the characteristic spatial volume of interaction 
and the power n is the number of initial state particles minus one. 
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In general, Eq. (0) generates a set of coupled nonlinear differential equations which 
can be solved numerically. However, instead of considering all the particle densities as 
unknowns, we can simplify the situation with the good approximation that the particle 
densities whose equilibrating chemical reaction rate is large compared to the Hubble 
expansion rate follow equilibrium densities of the form Eq. (f|). This, in fact, is the 
justification for our Boltzmann evolution's initial condition which is to start all species 
at equilibrium densities given by Eq. (|5|). With this expectation, we replace the X and 
the 7r densities in the Boltzmann equation with the equilibrium densities. We are then 
left to consider only the R° and the 7 densities as functions that require solutions. 

To understand which reactions will be most important in our system, we first recast 
Eq. (0) into the dimensionless form 

x_± Y (W jA ^ B] V^) U^ Ajl n x q ( UxeA Jt n Yp UkeB 3l Y k \ 

Yjdx j2 - H(x) { Y? Y i U keBji Y? ) [ } 

where Y T = n r /s, s is the entropy per comoving volume given by s ~ (27r 2 /45)g*m 3 /x 3 

(entropy conservation is assumed), and x = m/T. Note that we can interpret the 

numerator above H(x) to be the reaction rate per unit density of j's. For the purpose 

of illustration, suppose two reactions named a and b are governing the evolution of j 

particles and the reaction rates corresponding to them are labeled R a and Rb- The 

evolution equation in the form of Eq. (§) then becomes 

x dYj R a (x) , , R b (x) 

Y 3 -^ = ~ l^7 (ratlos a) - lw (ratlos b) (9) 

where the "ratios" refer to the terms consisting of density ratios. Suppose further that we 
are at a time when R a (x)/H(x) ^ 1 while R^{x)/ H{x) ^> 1. Then as long as the "ratios 
a" and "ratios b" are comparable in value, reaction a can be neglected during this period 
of evolution. Furthermore, if the final products of reaction b are in equilibrium, the j 
particle density will follow the equilibrium density as long as reaction b dominates. With 
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such reasoning, Ref. || argues that jR° <-> X and 77 <-> X reactions play a negligible 
role compared to <-» X, 771 <-> i?°7r, and 77? <-> i?° in keeping the R° and the 

7 densities in equilibrium near the time of 7 freeze out. In our present work, we shall 
neglect only the weakest of the relevant reactions, 7-R «-> X.Q 

The Boltzmann equations relevant to calculating the 7 abundance thus reduce to 
a pair of coupled differential equations containing terms corresponding to the set of 
reactions jnir <-> R°, jir <-» R°, 771 <-> i?°7r, R°R° <-> X, 77 <-> X: 



Yro (ix 





(10) 



i2w = (<W~^ F)<« + <W^^^ (12) 

The factor of X comes from summing over the isospins of 7T. In the next section, we 

argue that only 7r =l= should be included in 7V — > 7ri?°, resulting in X = 2.0 

Before we move on to discuss the transition rates, let us clarify the term "freeze out 

time" used in this paper, particularly in Section IV. In agreement with what will be 

revealed in the next section, suppose that the self-annihilation term in Eq. (|10"D can be 

neglected compared to the term associated with R tot . When R tot /H becomes much less 

than unity and continues decreasing sufficiently fast to keep the right hand side of Eq. 

(|T0|) much less than unity despite the increases in the magnitude of the term multiplying 

Rtot/H, the fractional change in Y~ becomes negligible. This is then a sufficient condition 

for the number of 7's becoming approximately constant (freezing out). We shall use the 

7 We have checked numerically that this reaction plays a negligible role. 
8 The choice N = 3 was implicit in the treatment in Ref. M. 
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term freeze out time to refer to the approximate time at which R to t/H becomes much 
less than unity. 

III. THE TRANSITION RATES 

Transition amplitudes for R°, 7, and pions depend on hadronic matrix elements of 
four-fermion effective operators of the form g'jqq, obtained by integrating out the squark 
degree of freedom. Since only a small number of fundamental short- distance operators 
underly all the transition amplitudes of interest, crossing symmetry can be used to relate 
transition amplitudes for some of the reactions. Due to the possibly strong momentum 
dependence of the amplitudes, however, this proves to be of limited utility. This is 
discussed in Section V. 

The particles R° and 7 are charge conjugation even and odd, respectively ||. Thus, 
if charge conjugation were a good symmetry of the interaction, pions coupling an R° 
to a 7 would have to be in a C-odd state. However, C invariance is violated by the 
mass-splitting between L and i?-chiral squarks (superpartners of the left and right chiral 
quarks). This mass splitting is a model dependent aspect of SUSY breaking. Fortunately, 
as we will see, our analysis is quite insensitive to the extent of C violation. 

We now present expressions for the transition rates to be used in Eq. fllO|) : (W-^^^V)^ 
(^o^ 2 }«, (W R0R0 ^ x V)n% , and (W~+ X V)r%. We shall see 
that the resulting expressions do not differ significantly from those of Ref. |§ even though 
the issue of charge conjugation symmetry is ignored in that reference. 

A. The 7 — R° conversion reaction 77? — > 7ri?° 
If charge conjugation invariance were exact, the neutral pion channel would be ab- 
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sent as it necessarily violates C. However, even if C is maximally violated, the condition 
cr(77r° — ► n°R°) <C cr(77r ± — > n^R ) still applies because the •yR°qq coupling is pro- 
portional to the quark charge, causing a first order cancelation to occur in the case of 
a neutral tt°. Thus we can ignore the neutral channel without serious impact on the 
quantitative results. In order to avoid any thermal averaging complications that may 
arise from threshold effects, we estimate the cross section for i?°7r ± — > 7r ±; y instead of 
the cross section for the inverse reaction. The cross section formula is the same as that 
given in ||: 

(W R0 ^ % ±~V) = (va ROw ) ~ 1.5 x l(T 10 r [/fe 4 C] mb. (13) 

The factor C contains the uncertainty due to possible resonance effects and hadronic 
physics. Ref. considered the range 1 < C < 10 3 . An analysis of the effect of the 
expected R n resonance shows that C can exceed this by an order of magnitude (see 
Appendix), but we shall not dwell on this since our conclusions are mostly insensitive 
to the exact value of any large enhancement. However, for C j \x\ < 1 the results are 
sensitive to the value of C j For reasons to be discussed in Section V, we also consider 
values of C as small as 1/20. 

Using Eq. <g) and nf « g j (m j T/(27v)) 3 / 2 e- m ^ T \ we then find 



7 " n~ 



= 9.17 x 10" 13 r 5/2 x~ 3/2 exp(-0.175 / u 8 1 a;)GeV 

xex P [-(r-l)x] WJ 2 ^s A C] . (14) 

Note in Eq. (|14D that parameters C and occur only together in the combination C / fig, 
which we take to lie in the range 

6.17 x 10" 4 < 10 4 (15) 



(4 
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in accordance with the limit on /is given by Table I and Eq. ([!]). 



B. The inverse decay reactions 

We now estimate the decay rate of R° and use Eq. ([?]) to obtain the inverse decay 
rate. If charge conjugation invariance is exact, two body decays of an R° to a 7 and 
a pseudoscalar meson (C=+l) are forbidden ||. In order to avoid reliance on a model 
of SUSY-breaking and its predictions for the extent of C-violation, we parameterize the 
branching fraction of an R° to 2- and 3-body final states by b 2 and 63, respectively. As 
in the 7 — R° conversion reaction, the neutral pion channel (R° — > 77r°7r°) can be safely 
ignored even if 63 ^> b 2 . When b 2 is not negligible, the two-body final states could be 77?° 
and 777. However, the matrix element-squared for R° «-»• 777 is about one-quarter of that 
for R° <->• 7T°7 0, and the rj final state is additionally suppressed by phase space. Hence 
we make an unimportant error by retaining only the two-body final state 77?°. Thus, the 
two reactions of interest are R° — > ; yir + ii~ and R° —>■ 777°, with branching fractions 63 
and b 2 , respectively. In this subsection, we show that our results depend only minimally 
upon the individual magnitudes of 62 and 63 because the Boltzmann equation depends 
only on the total decay width of the R° and b 2 + &3 ~ 1 (due to the relative phase space 
suppression of 4-body decays). 

The rates T(R° — > 777°) and T(R° —>■ 77r + 7r~) are obtained from the R° decay rate in 
Ref. @ by inserting 6 2 and 63 to get 

W#_^r- = r«o_^- = 2.0xlO- 14 ^(r)£(r-0.35(W 1 -l) GeV [^S 4 Bb 3 ] (16) 

and 

W R0 ^o = T R0 ^ = 2.0 x 1(T 14 F(r)6(r - 0.175/V - 1) GeV [^s 4 Bb 2 } (17) 

where J~(r) = r 5 (l — r -1 ) 6 , 9 is a step function employed to model the threshold of the 
decay channel, and the factor B reflects the overall uncertainty which we set to be in 
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the range 1/300 < B < 3. Having obtained the decay rate formula, we now use Eq. ([?[) 
to convert it to the inverse decay ratePI 



eq 

(^./W = ( r *p^-)^r ( 18 ) 

7 

= 2.0 x 10- u r 3 / 2 F(r)e- ir - 1)x 6(r - 0.350/i 8 1 - l)GeV[^ 5 4 56 3 ], (19) 
and similarly 



eq 

{W~ M0 V)n% = (T R o^ n o)^ (20) 

1 1 — 

7 

= 2.0 x 10- u r 3 / 2 F(r)e- {r - 1)x 6(r - 0.175/i 8 1 - l)GeV[/i^/i 5 4 56 2 ]. (21) 
Combining Eqs. (|P9|) and (|21"|) with Eq. (|12l), and using 5 2 + 6 3 « 1, we find 



= r io ^(6 2 , r, /i 8 ) + 2(W^ ± _ +7r±J? o^>< ,? ± (22) 

where T tot = 2.0 x lO-^r^^^e-^-^GeV^lfi^B} and 

[ 1 if r > 0.35^ 1 + 1 
g(b 2 ,r,ti 8 ) = < b 2 if 0.35/i 8 1 + 1 > r > 0.175/i 8 1 + 1 . 
[ otherwise 

The function g allows both the two and the three body decays when the R° is sufficiently 

heavy (r > 0.35/a 8 1 + 1) but forbids the three body channel when the R° mass drops 

below the two pion channel threshold. Thus, as long as the parameterization in Eqs. ([19]) 

and (|2lD is valid and the R° is massive enough (r > 0.35/i 8 1 + 1) to allow kinematically 

three body decays, our results are independent of 6 2 and 63, and hence the question of C 

A more accurate relationship between the thermal averaged decay rate and the non-thermal aver- 
aged decay rate is (T R0 _ t ~ 7r+ ^_ ) = T R0 _ i ~ n+n _Ki{rx) / K 2 {rx) where K v is the modified Bessel function 
irregular at the origin. Since the freeze out occurs typically between x — 20 and x — 30, the thermal 
averaged reaction rate will maximally deviate from the non-thermal averaged one when r = 1.1 and 
x = 20. In that case K\(jx) j K^irx) — 0.94 which is still an insignificant correction. Thus we neglect 
this complication in our calculations. 
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invariance. Therefore, considerations of C invariance is generally unimportant for large 
values of r. 

In these formulae, the squark mass parameter fi$ occurs only in combination with 
the uncertainty parameter B in the form B/fig. Using Table I and Eq. ([TJ), we limit the 
decay rate for a given //g and r to those values corresponding to the ranged] 

4.12 x 1(T 5 < — ^ 48- (24) 



C. Self- Annihilations and co-annihilations 
For the thermal averaged R° self-annihilation cross section, we use (va^ojio) = 31 A mb. 
This is extracted from the pp annihilation cross section in the comparable kinematic re- 
gion [[HJ with a factor A inserted to cover a possible difference between R°R° and pp 
annihilation, and to account for the uncertainty due to possible resonance enhancements 
and other hadronic effects. We take A to lie in the range 10~ 2 < A < 10 2 .Q Hence, the 

10 Because of our estimated upper and lower limit on each of the parameters C, B, and us separately 
(Eq. (|l5|), Eq. (|24|), and Table I), for a given value of -B/Ms' the allowed range of values for C j \i\ given 



by Eq. ( |15| ) must be supplemented with the condition 

(23) 





B 


Ms 


Cm ax 


B 


-E^max _ 


Ms 


-E$rnin _ 


Ms 



where from Eq. C max = 1000, C mm = 1/20, and B max = 3, and B mm = 1/300. 

11 Note the absence of the v 2 factor which appears in the familiar case of two identical Majorana 
spinors annihilating to a fcrmion-antifcrmion pair (e.g., 77 or gg — * qq). Like the 77 and gg states, 
the R°R° system must be antisymmetric by Fermi statistics, i.e., 1 Sq, Pi,.... However typical final 
states of R a R° annihilation (e.g., 3 pions) can have ^ quantum numbers, allowing s-wave annihilation. 
This is to be contrasted with the usual case that the final state is a fermion-antifermion pair. Since the 
sfermion-fermion-gaugino interaction conserves chirality, the ^ state in that case is helicity-suppressed 
and thus p-w&ve annihilation is necessary. This treatment departs from Ref. J9|, but does not lead to 
significantly different conclusions than the (w^o^o) = lOOAv 2 mb used there. 
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R° self-annihilation rate is given by 

(W B p B p_+ x V)n% = 10Ar 3/2 x~ 3/2 exp(-rx)[/i 3 ]GeV. (25) 

Note that although the R° self-annihilation rate is generally much larger than the other 
reaction rates before the 7 freeze out time, it is usually not strong enough to maintain 
R° in equilibrium abundance through the 7 freeze out time. This fact, not taken into 
account in Ref. ||, leads to differences in the results between that paper and the present 
analysis. 



The well known thermal average of the 7 self-annihilation cross section |l^Jl^JT6| JT 



can be approximated as |§ (vcr~) = 2.0 x 10 11 x 1 [/i§/ i 5] mD f° r our purposes, giving 
the transition rate 

<Wfi_ x V)v% = 3.3 x 10- 12 x- 5 / 2 exp(-x)[^^ 4 ]GeV. (26) 

Because the 7 self-annihilation becomes ineffective earlier than the R° self-annihilation, 
it contributes very little to our results. 

In summary, the reactions that will be important to our system of equations are 
R°R° <-> X, RPtt «-> 7vr, R° <-> 77T7T, and R° «-> 7tt. 



IV. GENERAL RESULTS 

In this section, we impose the cosmological constraint Q~h 2 < 1 on the integra- 
tion results of the Boltzmann equation to identify the allowed region of the parameter 
space and use the condition VL~h 2 > 0.01 to identify those parameters for which the 
photinos are significant dark matter candidates. The parameter space is spanned by 
r = M/m,B/i4,C/i4,lJifi, and A (see Table I and Eqs. (g), (0), (|IJ), and ([§)). For 
reasons of physical interest, we will present our results in terms of the 7 mass m and the 
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(ct( R° R° -X) v >= 31 mb 




0.6 0.7 0.8 



0.9 1 
m (GeV) 



1.1 1.2 1.3 1.4 



Fig. 1: For any given contour type, the left contour gives those values of M (R° mass) and the m 
(7 mass) for which Q,~h 2 =\ while the right one gives those for which Vt~h 2 = 0.01. The region above 
the left contour is ruled out by the present analysis. 

R° mass M instead of using r and /!§• We constrain the parameter space for the two 
extreme cases 62 = 1 and 62 = (maximal C violation and C conservation, respectively), 
but in general, the results are insensitive to the value of b 2 - As will be discussed below, 
among the parameters of the model, the relic abundance is most sensitive to the varia- 
tions of r. Using a maximum A of about 100, our analysis gives us an upper bound of 
r < 1.8. 

In Fig. 1, we show the Boltzmann equation integration results with exact C invariance 
(b 2 = 0). For any given contour type, the left contour represents the Q~h 2 = 1 contour 
while the right contour represents the VL~h 2 = 0.01 contour. The present analysis thus 
excludes the region above the left contour and constrains the masses to lie between a 
given contour type in order for the 7 to be a significant source of dark matter (defined 
by VL~h 2 > 0.01). In this figure, the parameter A multiplying the R° self-annihilation 
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Fig. 2: Same as Fig. 1 except for the BP self-annihilation cross sections. For the top figure, 
(va(B°B -> X)) = 3100 mb while for the bottom figure, (va{B°B a -> X)) = 0.31 mb. 
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cross section has been set to 1. 

Note that the values of r = M/m are insensitive to C ' j ' [i% > 1. This can be heuris- 
tically understood by the fact that as C / fx% increases, the freeze out time (the time 
at which the 7 — R° conversion reaction rate becomes negligible compared to H if the 
7 — R° conversion rate dominates over the inverse decay rate) approaches the time at 
which the R° self-annihilation rate becomes negligible compared to H. Thus, as Cj [i s 
increases, the photino abundance should approach the value for the limiting case when 
the R° self-annihilation rate becomes negligible before the freeze out time. When the 
R° self-annihilation rate becomes negligible, the number of SUSY particles are approxi- 
mately conserved-PI Thus, the 7 abundance is largely determined by the time at which 
R° self-annihilation becomes negligible in this limiting case. This time is determined by 
r and is independent of C//i|.0 

Because the extent of C violation affects the photino abundance only in the region 
0.14GeV < M — m < 0.28GeV, it is clear from Fig. 1 that only the long dashed contours 
(corresponding to small 7 — R° conversion and inverse decay rates) may depend on the 
extent of C violation. However, for this region, B / '//g is too small for the inverse decay 
reaction to play any significant role, and hence, our results are insensitive to the extent 
of C violation. Explicit numerical calculations confirm this. 

In Fig. 2, we show the effect of changing the magnitude of the R° self-annihilation 

cross section (by changing A in Eq. (p5p ). When we increase the magnitude from that 

of Fig. 1 by a factor of 100 (due to a possible resonance enhancement), the contours 

12 The 7 self-annihilation rate is already negligible by the the time the R° self-annihilation becomes 
negligible. 

13 This heuristic argument assumes that the R and 7 abundances approximately follow a function 
independent of C / fig until near the time that the SUSY particles become approximately conserved. 
This is of course true for the equilibrium abundance functions. 
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for C//i|, B/pLg > 1 shift leftwards, and when we decrease the magnitude by a factor of 
100, the same contours shift rightwards. In both cases, the contours corresponding to 
a small C j \x% and B / fig (corresponding to small 7 — R° conversion and inverse decay 
rates) remain essentially unchanged. This is expected since for the shifted contours, the 
inverse decay and the 7 — R° conversion reaction rates are large enough such that the 7 
abundance is sensitive to the time at which the R° self-annihilation becomes negligible 
(by the mechanism discussed before) while for the unchanging contours, the 7 abundance 
is determined nearly independently of the time that R° self-annihilation rate becomes 
negligible. When the inverse decay and the 7 — R° conversion reaction rates are very 
small as is the case for the unchanged contours, the 7 freeze out time that will lead to 
Q^h 2 > 0.01 is much earlier than the time when the R° self-annihilation reaction rates 
become negligible. Hence, near the 7 freeze out time, the R° self-annihilation reaction 
rate will dominate the dY R o/dx, the R° abundance will be nearly in equilibrium, and 
the dY~/dx will decouple from dY R o/dx, leading to a 7 freeze out value that is nearly 
independent of the R° self-annihilation reaction. Note also that when the time at which 
the R° self-annihilation becomes negligible is pushed away from the 7 freeze out time 
by increasing the R° self-annihilation cross section, the solid and the dotted contours 
become more sensitive to the value of C//z| as we expect from our heuristic discussion 
above. 

According to Fig. 2, the maximum value of r allowed by the condition that Vt~h 2 < 1 
is about 1.8. 

V. CROSSING RELATION 

The amplitudes determining the quantities (va) = (vcr R o w ±~ n ±) and T tot are related 
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through crossing symmetry if we associate T to t with the C conserving R° — > 77r + 7r~ tran- 
sition rate (i.e., if we set 62 = 0). To obtain a useful constraint from the crossing relation, 
and to implement the constraints following from the symmetries of the underlying the- 
ory, we derive in this section an approximate effective interaction Lagrangian. If the R n 
resonance is sufficiently far above threshold such that the R ^ — > 77r amplitude can 
be taken to be momentum independent for the purposes of the freeze-out calculation, a 
single parameter governs both (va R07V ±~ n ±) and T tot . This allows us to determine what 
ranges of R° lifetime are most favorable for cosmology, in the event the R^ is too far 
above threshold to have a significant impact. 

We first note that neglecting light quark masses as well as left-right squark-mass 
splitting in comparison to the squark masses, the four- Fermi effective operator governing 
R°7T <-> 77r can be written in the current-current form 

Hint = iK V Af 7 M A~ q 1 ^ + k a Af 7 M 7 5 A~ q 1 ^^ (27) 

where A~ and A~ are 4-component Majorana spinor fields for the gluino and photino, 
{a, i, j} are color indices, and the T a are 3x3 SU(3) matrices. This form follows 
because the underlying theory conserves the chirality of light quarks and their SUSY 
partners]^] allowing only current-current couplings for the quarks to appear. Approxi- 
mate degeneracy of the left-right squark masses then ensures parity conservation which, 
with Lorentz invariance, results in the form of Eq. (^). A direct calculation starting 



with the fundamental supersymmetric Lagrangian of course gives the form (|27|) and gives 



Ky = and ka = gse g e/ 'Ms 2 .[^| The vanishing of Ky is due to C-conservation, since the 

14 Chirality conservation of the light quarks and their squarks is an excellent approximation in all 
SUSY models proposed to date, for which left-right squark mixing is proportional to the mass of the 
corresponding quark. 

15 The strong coupling constant is denoted by gs and e q gives the electric charge of the quark in units 
of positron charge e. 
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term it multiplies is C-odd for Majorana fields A~ and A~. 

We are concerned with estimating matrix elements such as (i?°7r|7ii nt |77r). The most 
general form of the matrix element includes current- current terms, plus other terms 
which result from the fact that the R° is not pointlike and chirality flip can be induced 
by long-distance effects. However, since the R° is expected to be more compact than 
ordinary hadrons (as is observed for the ++ glueballQ, we neglect all but the current- 
current terms. Therefore we have C e s = ikRP^X-J^ where is a C-odd£] four- 
vector pion current determined by chiral perturbation theory, and k is of order 
The single-pion contribution to J M vanishes, and the two pion contribution is simply 
Jfi = i(ir dpTT — (9 M 7r)^7r). In general, k is a function of kinematic invariants, but far from 
resonances a constant should be a reasonable approximation. 

Using £ eS we can compute both (va) = (v<j R o n ±^~ w ±) and Y to t in terms of the single 
parameter k. Thus, for a given r and M, VL~h 2 is a function of the single parameter k. 
Likewise, values for {Q~h 2 , r, M} pick out a unique value of k, which in turn determines 
T to t- In Fig- 3, we assume Q~h 2 = 0.25 (cosmologically "favored" value) and give the 
R° lifetime for a range of r and R° mass. We stress that these results are only indica- 
tive of the actual lifetime-mass-relic density relation, since the most general effective 
Lagrangian depends on additional parameters which we neglect here. Furthermore if the 
R n resonance is sufficiently close to threshold to produce an enhancement effect, there 
is no simple relation among Q~h 2 , r, and the R° lifetime, and Fig. 3 is not relevant. 
It is encouraging that the R° lifet imes required to give the "correct" relic density is 
compatible with predictions and also compatible with experimental limits ||. 

It is also of interest to extract the values of B and C implied by VL~h 2 = 0.25; this 
is shown in Fig. 4. Overall, these results suggest that the inverse decay reaction may 
16 D. Weingarten, private communication. 

17 Because the R° and 7 have opposite C quantum numbers ||. 
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Fig. 3: The R° lifetime that implies a cosmological photino abundance of £!~/i 2 = 0.25 is plotted 
as a function of the R° mass and its ratio r to the photino mass. A model Lagrangian has been used to 
determine the crossing relation between the 7 — R° conversion amplitude and the R° decay amplitude. 




Fig. 4: The -B/Ms an d C/^s values corresponding to the contours shown in Fig. 3 are plotted. The 
typical suppression of Cj \i% with respect to B j [i A s reflects the fact that both the 7 — R° conversion and 
the inverse decay reactions have comparable rates in our simple model which does not take into account 
possible resonance enhancements. 91 



not be entirely negligible in determining the photino abundance if there is no resonance 
enhancement in the scattering reaction. 

VII. SUMMARY AND CONCLUSIONS 

In this article, we have investigated the cosmological constraints on the physics of light 
photinos and gluinos. A full treatment of the Boltzmann equations governing the photino 
freeze out has been carried out, considering the total R° width and R°tt — > 7V scattering 
cross section as independent quantities. We find that to avoid photino abundances 
inconsistent with cosmology, the ratio r of R° mass to 7 mass must be less than about 
1.8. We checked that if the R° is the LSP, its annihilation is too efficient for it to account 
for the observed dark matter density. 

We also developed an approximate effective Lagrangian description of the i?°7r <-» 
77r amplitude, neglecting possible C-violating and chirality-violating effects. If the R v 
resonance is far above threshold, C e s is specified by a single parameter governing both 
the total R° width and R°tt — > 77? scattering cross section. Assuming that the universe 
is at its critical density with photinos constituting most of the dark matter fixes this 
parameter for given R° and 7 masses. We therefore obtain the cosmologically favored 
lifetime range of the R° as a function of its mass (shown in Fig. 3) in the absence of a low 
lying R n resonance. The lifetime will be increased compared to the values given in Fig. 
3 when the R n resonance enhances the cosmological importance of the scattering cross 
section in comparison to the inverse decay. Although the limitations in this estimate 
must not be forgotten, it is encouraging that the range thus determined, r > 10 _10 s, 
is compatible with experimental limits [|J. Much of this range of lifetimes should be 
accessible to direct observation in upcoming experiments 0. 
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In closing, we note that detectability of relic dark matter is different for light 7's than 
in the conventional heavy WIMP scenario for two reasons. Firstly, the usual relation 
between the relic density and the WIMP-matter scattering cross section only applies 
when the relic density is determined by the WIMP self-annihilation cross section, whereas 
in the light photino scenario it is determined by the 7 — R° conversion cross section, 
R° self-annihilation cross section, and the density at freeze out. Secondly, WIMP 
detectors have generally been optimized to maximize the recoil energy for a WIMP mass 



of order 10 to 100 GeV. Goodman and Witten in Ref. [18] discuss 7 detection through 



7-nucleon elastic scattering. Using Eq. (3) of Ref. ]T8| and the parameters discussed 
here, one finds that event rates range somewhere between 10 -3 and 10 events/(kg day). 
Unfortunately even if the event rate were larger, observation of relic light photinos would 
be difficult with existing detectors because the sensitivity of a generic detector is poor 
for the less than 1 GeV mass relevant in this case. 
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APPENDIX 



In this Appendix, we give the formalism for treating the resonance enhancement of 
the R°iT — > 77? cross section, using a Breit-Wigner form for the resonance. This permits 
us to assess the plausibility of the original range used in Ref. 0, C < 1000. We find that 
the effective value of C could be significantly larger than the originally estimated upper 
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bound, but this is only relevant if (va R o R o) is large enough that R 0, s remain in thermal 
equilibrium until after photino freeze out. As discussed in Section IV, this is not the case 
for large C, given our estimated (va R o R o). However if there were a h glueball near 
R°R° threshold, the R 0, s could stay in equilibrium to a lower temperature, and make it 
necessary to include resonance effects for both self-annihilation and 7 — R° conversion 
processes. We treat below the modeling of a resonance in the R n — > 77? reaction; the 
treatment of a resonance in the R°R° self-annihilation cross section is a close parallel. 

The resonance relevant to the R°tt — ► 77r reaction is called R^ which is composed at 
the valence level of g, qi, and q 2 (where q^s are u and d quarks). To study the maximum 
enhancement, we consider the R n mass to be close to the R° mass. We also consider here 
only the charged R n 's since we are concerned with charged pion scattering (see Section 
III). Furthermore, because the s-wave contribution dominates, we restrict ourselves to 
the J = 1/2 state. 

We write the resonant contribution to the R°tt — > 77? cross section as 



where p cm is the center of mass three-momentum of the incoming particles, s is the square 

of center of mass energy, is the mass of R n , the T(A — > BC)'s are momentum (s) 

dependent widths to the incoming and outgoing channels^ and T to t is the momentum 

dependent total width of the R n . Thus, Eq. fll3|) becomes (va R o n ) = (vcr nonres )-|-(f(Tres) 

where (v anonres) is the formula given in Eq. (p~3[ ) with C set to a value of order 1. Since 

we are concerned with the maximum cross section resulting from the resonance, we 

are focusing on the region of parameters for which the non-resonant cross section is 

unimportant (i.e. (wnonres) < (^crres))- 

18 Note that C poses no relevant constraints on the decays of the charged so both R w ± — > R ^ 
and R^± — ► jir^ are allowed even though the R° and 7 have opposite C eigenvalues. 



Ores — 



Pcm 



An m 



2 R r(R w -> 77r)r(K ^ R°k) 



(28) 
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The kinematic momentum dependence of T(A —>■ BC) can be seen by expressing it 
in terms of a solid angle integral over the invariant amplitude squared |.M| 2 : 

T(A - BC) = [p cm (s)/(32vr 2 m fcV ^)] j dQ\M(A - BC)\ 2 . (29) 



Here p cm (s) = l/2y (s — (tub + rnc) 2 )( s — ( m B — mc) 2 )/ 's is the center of mass frame 
three-momentum of the decay products B and C. Defining Att^b = I dQ\Ai(A —>■ BC)\ 2 
and assuming that T tot ~ T(R n — > R°n), the three independently adjustable parameters 
for the resonance are taken to be m Rn ,^ R o, and £~. 

Since R n — > R° + 7r is a strong decay, we can take its matrix element to be similar 
to the matrix element for some known strongly decaying resonance whose decay has 
no angular momentum barrier, for instance the /o(1370) whose total width is 300-500 
MeV|5|. Thus we use £ R o « 16tt r(/ (1370)) m/ (i 3 7o) = 7.4 GeV 2 . 

To determine £~ we estimate the ratio of R w — > 77? and R w — > i?°7r matrix elements 
by keeping track of the factors entering the short distance operator responsible for R° — > 



7T7T7, namely 7gg<2'. We use the i?° mass, M, to set the scale. This gives 

= A e2 ^T^ - ^ 4 >< 10 " 10 ^ ^8 % 4 - (30) 



We define C e to be the effective value of C in Eq. (|T3|) that would reproduce the f2~/i 2 
calculated using the present resonance model for a given set of resonance parameters, 
and taking A large enough to keep the R° in equilibrium abundance until after photino- 
freeze out. To calculate the thermal average (wres) of the resonant cross section, we use 
a non-relativistic approximation which is within a factor of two or better of the exact 
average. It can be expressed in terms of the one- dimensional integral 



(vares) = + J) f*( S («)^ (31) 

where s((3) « 0.64/i 2 (r + 0. 175 //i 8 )(r + 0.175//x 8 + 2/5) GeV 2 . 
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Fig. 5: The effective resonance enhancement factor C e /fj,g is shown as a function of r, m, and m^. 

C e is largest when the resonance is near the threshold of the R°ir channel, because near 
the threshold T(R n — > R°ti) is phase space suppressed in comparison to T(R n — > 771-) 
and the peak value of the Breit-Wigner cross section is proportional to ~ r(i? 7r — > 
77r)/r(_R 7r — > i?°7r). However, because the width of the resonance vanishes as mR n 
approaches threshold, the thermal average integral of the Breit-Wigner cross section 
does not grow arbitrarily large. 

In order to assess the plausibility of the original range used in Ref. 0, we plot (Fig. 
5) C e //x| as a function of r, m, and with A = 1 (in Eq. (|30|)) and /is = 1/2. In 
all of the m and r cases shown, C e //i| > 2 x 10 5 (or equivalently C e > 10 4 ) when 
— M — m,r < 70 MeV. Since the mass splitting can easily be less than 70 MeV, we 
see that the C range used in Ref. would be inadequate, were self-annihilation to be 
significantly larger than the non-resonant estimate adopted here. 
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